Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

Prediction and functional analysis of native disorder in proteins from the Three Kingdoms of Life

An automatic method for recognizing natively disordered regions from amino acid sequence is described and benchmarked against predictors that were assessed at the latest critical assessment of techniques for protein structure prediction (CASP) experiment. The method attains a Wilcoxon score of 90.0, which represents a statistically significant improvement on the methods evaluated on the same targets at CASP. The classifier, DISOPRED2, was used to estimate the frequency of native disorder in several representative genomes from the three kingdoms of life. Putative, long (>30 residue) disordered segments are found to occur in 2.0% of archaean, 4.2% of eubacterial and 33.0% of eukaryotic proteins. The function of proteins with long predicted regions of disorder was investigated using the gene ontology annotations supplied with the Saccharomyces genome database. The analysis of the yeast proteome suggests that proteins containing disorder are often located in the cell nucleus and are involved in the regulation of transcription and cell signalling. The results also indicate that native disorder is associated with the molecular functions of kinase activity and nucleic acid binding.

Automatic prediction of functional site regions in low-resolution protein structures

World-wide structural genomics initiatives are rapidly accumulating structures for which limited functional information is available. Additionally, state-of-the art structural prediction programs are now capable of generating at least low resolution structural models of target proteins. Accurate detection and classification of functional sites within both solved and modelled protein structures therefore represents an important challenge. We present a fully automatic site detection method, FuncSite, that uses neural network classifiers to predict the location and type of functionally important sites in protein structures. The method is designed primarily to require only backbone residue positions without the need for specific side-chain atoms to be present. In order to highlight effective site detection in low resolution structural models FuncSite was used to screen model proteins generated using mGenTHREADER on a set of newly released structures. We found effective metal site detection even for moderate quality protein models illustrating the robustness of the method.

Functional benefits of predator species diversity depend on prey identity

1. Determining the functional significance of species diversity in natural enemy assemblages is a key step towards prediction of the likely impact of biodiversity loss on natural pest control processes. While the biological control literature contains examples in which increased natural enemy diversity hinders pest control, other studies have highlighted mechanisms where pest suppression is promoted by increased enemy diversity. 2. This study aimed to test whether increased predator species diversity results in higher rates of predation on two key, but contrasting, insect pest species commonly found in the rice ecosystems of south-east Asia. 3. Glasshouse experiments were undertaken in which four life stages of a planthopper (Nilaparvata lugens) and a moth (Marasmia patnalis) were caged with single or three-species combinations of generalist predators. 4. Generally, predation rates of the three-species assemblages exceeded expectation when attacking M. patnalis, but not when attacking N. lugens. In addition, a positive effect of increased predator species richness on overall predation rate was found with M. patnalis but not with N. lugens. 5. The results are consistent with theoretical predictions that morphological and behavioural differentiation among prey life stages promotes functional complementarity among predator species. This indicates that emergent species diversity effects in natural enemy assemblages are context dependent; they depend not only on the characteristics of the predators species, but on the identity of the species on which they prey.

Predicting functional gene links from phylogenetic-statistical analyses of whole genomes

An important element of the developing field of proteomics is to understand protein-protein interactions and other functional links amongst genes. Across-species correlation methods for detecting functional links work on the premise that functionally linked proteins will tend to show a common pattern of presence and absence across a range of genomes. We describe a maximum likelihood statistical model for predicting functional gene linkages. The method detects independent instances of the correlated gain or loss of pairs of proteins on phylogenetic trees, reducing the high rates of false positives observed in conventional across-species methods that do not explicitly incorporate a phylogeny. We show, in a dataset of 10,551 protein pairs, that the phylogenetic method improves by up to 35% on across-species analyses at identifying known functionally linked proteins. The method shows that protein pairs with at least two to three correlated events of gain or loss are almost certainly functionally linked. Contingent evolution, in which one gene's presence or absence depends upon the presence of another, can also be detected phylogenetically, and may identify genes whose functional significance depends upon its interaction with other genes. Incorporating phylogenetic information improves the prediction of functional linkages. The improvement derives from having a lower rate of false positives and from detecting trends that across-species analyses miss. Phylogenetic methods can easily be incorporated into the screening of large-scale bioinformatics datasets to identify sets of protein links and to characterise gene networks.

Anticipatory activation in the amygdala and anterior cingulate in generalized anxiety disorder and prediction of treatment response

OBJECTIVE: The anticipation of adverse outcomes, or worry, is a cardinal symptom of generalized anxiety disorder. Prior work with healthy subjects has shown that anticipating aversive events recruits a network of brain regions, including the amygdala and anterior cingulate cortex. This study tested whether patients with generalized anxiety disorder have alterations in anticipatory amygdala function and whether anticipatory activity in the anterior cingulate cortex predicts treatment response. METHOD: Functional magnetic resonance imaging (fMRI) was employed with 14 generalized anxiety disorder patients and 12 healthy comparison subjects matched for age, sex, and education. The event-related fMRI paradigm was composed of one warning cue that preceded aversive pictures and a second cue that preceded neutral pictures. Following the fMRI session, patients received 8 weeks of treatment with extended-release venlafaxine. RESULTS: Patients with generalized anxiety disorder showed greater anticipatory activity than healthy comparison subjects in the bilateral dorsal amygdala preceding both aversive and neutral pictures. Building on prior reports of pretreatment anterior cingulate cortex activity predicting treatment response, anticipatory activity in that area was associated with clinical outcome 8 weeks later following treatment with venlafaxine. Higher levels of pretreatment anterior cingulate cortex activity in anticipation of both aversive and neutral pictures were associated with greater reductions in anxiety and worry symptoms. CONCLUSIONS: These findings of heightened and indiscriminate amygdala responses to anticipatory signals in generalized anxiety disorder and of anterior cingulate cortex associations with treatment response provide neurobiological support for the role of anticipatory processes in the pathophysiology of generalized anxiety disorder.

Proteogenomics and in silico structural and functional annotation of the barley powdery mildew Blumeria graminis f. sp. hordei

Blumeria graminis is an economically important obligate plant-pathogenic fungus, whose entire genome was recently sequenced and manually annotated using ab initio in silico predictions [7]. Employing large scale proteogenomic analysis we are now able to verify independently the existence of proteins predicted by 24% of open reading frame models. We compared the haustoria and sporulating hyphae proteomes and identified 71 proteins exclusively in haustoria, the feeding and effector-delivery organs of the pathogen. These proteins are ‘significantly smaller than the rest of the protein pool and predicted to be secreted. Most do not share any similarities with Swiss–Prot or Trembl entries nor possess any identifiable Pfam domains. We used a novel automated prediction pipeline to model the 3D structures of the proteins, identify putative ligand binding sites and predict regions of intrinsic disorder. This revealed that the protein set found exclusively in haustoria is significantly less disordered than the rest of the identified Blumeria proteins or random (and representative) protein sets generated from the yeast proteome. For most of the haustorial proteins with unknown functions no good templates could be found, from which to generate high quality models. Thus, these unknown proteins present potentially new protein folds that can be specific to the interaction of the pathogen with its host.

Proteins and their interacting partners: an introduction to protein–ligand binding site prediction methods

Elucidating the biological and biochemical roles of proteins, and subsequently determining their interacting partners, can be difficult and time consuming using in vitro and/or in vivo methods, and consequently the majority of newly sequenced proteins will have unknown structures and functions. However, in silico methods for predicting protein–ligand binding sites and protein biochemical functions offer an alternative practical solution. The characterisation of protein–ligand binding sites is essential for investigating new functional roles, which can impact the major biological research spheres of health, food, and energy security. In this review we discuss the role in silico methods play in 3D modelling of protein–ligand binding sites, along with their role in predicting biochemical functionality. In addition, we describe in detail some of the key alternative in silico prediction approaches that are available, as well as discussing the Critical Assessment of Techniques for Protein Structure Prediction (CASP) and the Continuous Automated Model EvaluatiOn (CAMEO) projects, and their impact on developments in the field. Furthermore, we discuss the importance of protein function prediction methods for tackling 21st century problems.